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Abstract 

The coupling constants go* Dn and gs*Bn °f the- D*Dir and B*Bir strong interac- 
tions are computed in the effective quark model with chiral symmetry incorporating 
Heavy quark effective theory (HQET) and Chiral perturbation theory at the quark 
level (CHPT)g with linear realization of chiral U(3) x U(3) symmetry. We predict 
9d*d-k = 5.2 and g§*^ n = 14.3. The 1/M C corrections to the constant gD*Dn are 
calculated and make up less than 1%. The contribution of 1/Afj, corrections to the 
constant g^'B-K ls much less. 
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1 Introduction 



The investigation of strong interactions D*Dn and B*Bir runs parallel to the investigation 
of a dominance of vector D* and B* mesons in the semileptonic decays D — > 7r£i>£ and B — > 
7i£i>e related to the strong transitions D — > 7r and — > 7r, respectively, the form factors of 
which can be measured experimentally. In turn, the computation of the coupling constant 
9d* Dn °f the D*D7r interaction has a self-determined meaning, since the transitions D* — > 
Dn do exist on-mass-shell of interacting mesons. The probabilities of the decays B(D* + — > 



D°ir+) = (68.1 ± 1.3)%, B(D* + -> D+vr ) = (30.8 ± 0.8)% and B(D*° -> flV) = 



(63.6 ± 2.8) % have been measured recently by CLEO Collaboration [1]. 

The numerical values of Qb*d-k and g B * B7I are distributed from gD*Dw — 7 ±2, Qb-Stt = 



15±4 [2] and(? D * D7r = 6.3±1.9,^b w = 14±4 [3] up to = 12.5±1.0, 0b.b w = 29±3 



[4]. For the more complete list of references to the numerical values of the coupling 
constants go*Dn and gB*m we relegate readers to Table 1 of Ref.[4]. Most of these results 
have been obtained within QCD sum rules approach [5]. 

Recently [6] in the effective quark model with chiral U(3) x U(3) symmetry incor- 
porating Heavy quark effective theory (HQET) [7-9] and Chiral perturbation theory at 
the quark level (CHPT) 9 with linear realization of chiral U(3) x U(3) symmetry [10] the 
coupling constant gD*Dir has been computed at leading order in large c-quark mass M c 
expansion with the following value gD*Dn = gD*+D°Tr+ — 5.2 ± A. Here A is a theoretical 
uncertainty of our effective quark model which on the estimate of [6] makes up about 50%. 
However, for all cases of the application of HQET, supplemented by (CHPT) g , to the de- 
scription of strong interactions of heavy-light mesons [11-16] the resultant agreements 
between theoretical and experimental data turn out much better. 

A possible substantial change of the gD*Dw coupling constant value in HQET, sup- 
plemented by (CHPT) g , one can expect at next-to-leading order in large M c expansion. 
Indeed, as has been shown in Ref.[15] such correcitions to the form factors of semileptonic 
.D-meson decays can make up more than 30%. 

This letter is to compute qd*d-k in HQET, supplemented by (CHPT) g , to next-to- 
leading order in large M c expansion and to extend the obtained results on the computation 
of the coupling constant gB*B-w defining the B*B7r strong interaction. 

However, we have found that the contribution of 1/M C corrections to gD*Dn makes up 
less than 1%. In turn 1/M b corrections, where M b is the mass of the 6-quark, to g b* b-k is 
much less. The numerical value oi g B * B ^ calculated at leading order in large M b expansion 
equals: gs*B-K — 14.3. Thus, our results turn out to be compared reasonably well with 
the predictions of Refs.[2,3]. Therefore, in complete agreement with Refs.[2,3] we predict 
much less constibutions of vector mesons D* and B* to the form factors of semileptonic 
decays D — > n£u£ and B — > ixlvt than it can be expected due to Ref.[4]. 



2 Calculation of coupling constants at leading order 
in heavy quark mass expansion 



The coupling constant gD*Dw = gD*+D"w+ of the D* 
matrix element 



D°tt + decay is defined by the 



< out;ir + (q)D°(p)\D*(Q);in >= 



(2.1) 
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with 

M(D*+(Q) -> D°(p)7r+(q)) = g D * Dn q ■ e(Q), (2.2) 

where i£j (i = tt + ,D° or D* + ) is the energy of the i-meson, and V is the normalization 
volume, then e u (Q) is the 4-vector of the D* + -meson polarization. Note that by definition 
the coupling constant go*D-K in Eq. (|2.2j) involves a factor 2 with respect to that defined 
in Ref.fl] (see Eq.(8) of Ref.[6]). 

The computation of the matrix element M(D* + (Q) — > D°(p)-n- + (q)) at leading order 
in large N, number of quark colours, and M c expansion has been carried out in Ref. [6]. 
The result reads 

M(D* + (Q) -> D°(p)ir + (q)) = V2g ^ MD _ MD * £n(^ q . e(Q). (2.3) 
This gives one the coupling constant 

g D *+D^ = gD*Dir = in\ — I = 5.2, (2.4) 



where we have set M c = M D = 1.86 GeV and M D . = 2.00 GeV [1,2], and g n + = Att/VN 
[6,11-16] at iV = 3, then v' = 4A = 2.66 GeV [6,11-16] and A is the cut-off in Euclidean 3- 
momentum space connected with A x , the scale of spontaneous breaking of chiral symmetry 
(SBxS) in (CHPT) 9 , via the relation A = K x /V2 = 0.67 GeV at A x = 0.94 GeV [10]. 
The cut-off A appears due to the computation of heavy-light constituent quark loops 
describing strong low-energy interactions in HQET and (CHPT) g at scales p < A x . 

The coupling constant g^s-ir can be obtained from Eq.(|2.4|) by the change M c — > M b , 
M D -> M B and M D * -> M B * 

g B *B« = zr. tn[ — = 14.3, (2.5) 



TV v' \4m J 

where we have set M B « = M B = M b = 5.3 GeV [6]. 

Our predictions for the coupling constants go*Dn = 5.2 and g B * B7T = 14.3 agree with 
those given in Refs.[2,3]. 



3 Coupling constants to next— to— leading order in heavy 
quark mass expansion 

Now we should proceed to the analysis of next-to-leading order corrections in heavy 
quark mass expansion. We suggest to start with the computation of the 1/M C corrections 
to the coupling constant gD*Dn- In our approach, where the amplitudes are described in 
terms of constituent quark loops with virtual momenta restricted by the SB^S scale A x 
M c , 1/M C corrections can be computed by direct expansion of the c-quark Green function 
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in the momentum representation [15,16]. Since the virtual momentum k is restricted by 
the SB^S scale A x <C M c , we can expand the c-quark Green function in powers of k/M c 
and, holding only the terms of order 0(1/M C ), obtain [15,16] 



1 



1 +v 



M c -k 
1 

~2M r 



p 



v ■ k + iO 



1 + v 



k' 2 



v ■ k + i 



(v-k + iOf 



o\ 



M 2 



(3.1) 



where, as usually, we have set 
D°(p)iT + (q)) can be given by 



M D v^ 



M c v». 



Thus, the amplitude M(D* + (Q) 



M(D*+(Q) - D°(p)ir + (q)) = 

= M (0) (/r+(Q) -> D°(p)ir + (q)) + M W (D* + (Q) -> £>°(p) (3-2) 
7r + (g)), (3.3) 

where the term M^\D* + (Q) -> D°(p)7r+(g)) is given by Eq. (^3|) , while MW(L>* + (Q) -»■ 
D°(p)7r + (g)) reads 



1 



tr< ; :7 ^7 

[m — k + q m — k 



v ■ k + i 



y/2 2M C 16vr 



I 7m 



(v-k + iOf 



(3.4) 



After the calculation of the trace over Dirac matrices we reduce the r.h.s. of Eqj 
the form 



to 



M^{D* + {Q)^D\ P )-K + {q)) 



3d 9d* 



9^ 



N 



v/2 2M C 16tt 2 
d 4 k 4(m 2 - k 2 ) k fl + 2k 2 q fl 1 



e"(Q) x 



7i 2 i \m 2 



d 4 k 



(k — q) 2 — i 0][m 2 
2 m k 2 q^ 



k 2 - z 0] [v • fc + z 0] 
1 



7r 2 i [m 2 — (k — q) 2 — i 0] [m 2 



fc 2 -zO] [u • & + z'0p 



(3.5) 



Following the prescription of applied in our effective quark model [6,10-16] (see also Ref. 
[17]) for the calculation of constituent quark loop diagrams we have held only the leading 
divergent contributions dropping out the contributions independent of the cut-off. This 
prescription is nothing more than a trivial phenomenological description of quark confine- 
ment. Since the appearance of free quarks in the intermedite states of low-energy hadron 
interactions described by the constituent quark loop diagrams with constant constituent 
quark masses is related to the imaginary parts of these diagrams, the trivial prohibition 
of such an appearance is to keeping the divergent contributions of quark loop diagrams 
related to the real parts of quark diagrams. Moreover, the divergent parts of quark dia- 
grams, being polynomials in powers of the momenta of interacting hadrons, provide local 
effective low-energy interactions of hadrons in agreement with effective Chiral Lagrangian 
approach. In terms of the cut-off and the constituent quark mass this procedure restores 
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fully all effective coupling constants of low-energy interactions of low-lying mesons rep- 
resented at the hadron level by effective Chiral Lagrangians [10]. 

The integral containing k^ in the integrand has the following Lorentz structure: A(v ■ 
Q.1 Q 2 ) v ti + B{y ■ q, q 2 ) q^. The contribution of the term linear in is of order 0(1/M 2 ) 
due to Q,j,e^{Q) = and should be dropped out. Therefore, only the terms linear in 
should be held and give the contributions leading in chiral expansion. The integration 
over k gives one 



/ 



d 4 k 



4(m 2 - k 2 )k^ + 2k 2 q^ 



1 



TT 2 i [m 2 — (k — q) 2 — iO] [m 2 — k 2 — iO] [v ■ k + iO] 

d 4 k 2mk 2 q ll 1 

7i 2 i [m 2 — (k — q) 2 — iO] [m 2 — k 2 — iO] [v ■ k + iO] 2 

( r> \ 
= —8 mm — 
\4m J 



(3.6) 



For the computation of the intergals in Eq. 
d 4 k 
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we have used the integrals 



— V 



4£n\ 



'f) 

4m J 



(3.7) 



We should accentuate that no logarithmically divergent contributions have been omitted 
in the linearly divergent momentum integrals. Indeed, for example, one can show that 



d 4 k 



1 



Tc 2 i [m 2 — k 2 — i 0] [v ■ k + % 0] 

. „ A 
—4 A + 4 m arctg ► — v + 2%m 
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k 2 + m 2 
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Thus, following our phenomenological description of confinement, demanding to drop out 
the terms independent of the cut-off, and collecting all divergent terms of order 0(1/M C ) 
we obtain 



M«(£>*+(Q) - D°{p)TT + {q)) 



V2 



9-K+ 



^M D M D , 



Am ) 12M C 



q-e(Q). 



(3.9) 



This gives the coupling constant gD*Dn to next-to-leading order in large M c expansion 

9D*Dtt 



V2 



9t 



a/ Mo Mr,* 



v' 



5.2 x (1 + 7.3 x 10 



-3\ 



v 4m J [M c \Am) 



\2M r 



5.24. 



(3.10) 



Thus, we have found that 1/M C corrections to the constant go*Dn make up less than 
1% and can be neglected. The accounting of the logarithmically divergent constributions 
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together with the linearly ones does not go beyond the scape of the approach. Indeed, for 
the correct calculation of the kinetic terms of the low-lying scalar mesons related to the 
kinetic terms of the low-lying pseudoscalar ones via chiral transformations one needs to 
take into account both quadratically and logarithmically divergent constributions [10]. 

For the computation of the coupling constant gB*Bw we ) following the procedure ex- 
pounded above, get 



9B* Bit 



V2g^ 
14.3, 



in 



v' \ m ( v' \ 
Am J M b \Am J 



12M h 



(3.11) 



where the numerical value is obtained at M B , = M§ = Mf, = 5.3 GeV [1]. Thus, we 
predict negligibly small next-to-leading corrections in heavy Mb expansion to the coupling 
constant gs*B-K — 14.3. 



4 Conclusion 

When applying the effective quark model with chiral U(3) x U(3) symmetry incorporating 
HQEF and (CHPT) g to the computation the coupling constants Qd'D-k and g^Bn to next- 
to-leading order in heavy-quark mass expansion, we have computed at leading order in 
large N expansion and in the chiral limit the numerical values go'Dw = 5.2 and g&Bn = 
14.3 compared reasonably well with the results obtained by Colangelo et al. [2] within 
QCD sum rules approach and by Dosch and Narison [3] within full QCD. Then, we 
have shown that next-to-leading corrections in large heavy-quark mass expansion give 
negligible contributions and the main values of the coupling constants are induced at 
leading order. Our results as well as those predicted by Colangelo et al. [2] and Dosch 
and Narison [3] are a factor of 2 less compared with those obtained by Belyaev et al.[4]. 
Therefore, in comparison with Ref. [4] we predict much less contribution of vector mesons 
D* and B* to the form factors of the semileptonic decays D — > ixtvt and B — > ixtvi. 
The former is very important for the understanding of a mechanism of the semileptonic 
decays D — > -nivji and B — > tt£u£. It should be stressed that our numerical results contain 
a theoretical uncertainty ±A, the rough estimate of which gives A ~ 50% [6]. However 
for all apllications of HQET and (CHPT)^ [6,11-16] the resultant agreements between 
theoretical and experimental data turn out much better. 
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